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ABSTRACT 

We present the numerical simulations for an electron-beam-driven and loss-cone-driven electron- 
cyclotron maser (ECM) with different plasma parameters and different magnetic field strengths for a 
relatively small region and short time-scale in an attempt to interpret the recent discovered intense 
radio emission from ultracool dwarfs. We find that a large amount of electromagnetic field energy can 
be effectively released from the beam-driven ECM, which rapidly heats the surrounding plasma. A 
rapidly developed high-energy tail of electrons in velocity space (resulting from the heating process of 
the ECM) may produce the radio continuum depending on the initial strength of the external magnetic 
field and the electron beam current. Both significant linear polarization and circular polarization of 
electromagnetic waves can be obtained from the simulations. The spectral energy distributions of the 
simulated radio waves show that harmonics may appear from 10 to 70^ pc (^ po is the electron plasma 
frequency) in the non-relativistic case and from 10 to 600^0 in the relativistic case, which makes 
it difficult to find the fundamental cyclotron frequency in the observed radio frequencies. A wide 
frequency band should therefore be covered by future radio observations. 

Subject headings: magnetic fields - radio continuum: stars - stars: low-mass, brown dwarfs - polariza- 
tion - masers 



1. INTRODUCTION 

Ultracool dwarfs (UCDs) are those objects with 
spectral type later than M7 and low luminosity. Re- 
cent observations of 193 UCDs reveal that 12 UCDs 
produce intense radio emi ssion with flux den si ties up 
to hu ndr eds of /xJy (e .g. IBerger et al.l (l200ll): IBerger 



2002); IBerger et al l 



Burgasscr & Putman 



20051): IHallinan et alj (|2006D : lAntonova et al.l (|2008() 



IMcLean et all t)2011l )V Some of them show a highly 
circularly polarized radio pulse with regular periods , 
and a flux d ensity up to 15 mJy (|Hallinan et alj [20071 : 
IBerger et all[2009] ). 

In a range of past studies, these radio features of 
UCDs have been presented as a function of magnetic 
field, spectral type, rotation, age, binarity . and associ- 
ation wi th the X-ray and Ha emission (e.g. IBerger et afl 
(pOOl : IHallinan et all (|2006t) : lAntonova et all (|2008l) : 



IBerger et al] (|2010t) : IMcLean et all (f201lT > ). The regu- 
lar periods of radio pulses from TVLM 513-46546 and 
the L dwarf binary 2MASSW J0746425+200032 indicate 
that the radio activit y of UCDs is strongly i n conjunction 
with their rotation (|Hallinan et al.l 120071 : IBerger et alj 
2009) which is one of the crucial factors to influence the 
magnet ic field by a d ifferential rotation induced dynamo 
theory ([Parker! Il955h although oth er mechanisms, such 
as a turbulence-induced dynamo (|Durnev et al.1 11993D 
(for s mall-scale fields) or a 2 dynamo ([Chabrier fc K iikcr 
2006) (for large-scale fields), can also contribute to the 
magnetic field. The topology of the magnetic field on 
UCDs may be understood as a dipole due to the narrow 
bunching of multiple pulses of bot h left- and righ t- 100% 
polarization ([Hallinan et al.ll2007t) . IBerger et all ([2009D . 
however, has suggested that the field topology maybe 
more complex - due to a 1/4 phase lag of the radio pulses 
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compared to Ha. This is based on the assumption that 
the emission is parallel to the magnetic field, but if the 
emission is perpendicular to the field, then the 1 /4 phase 
lag is in agreement with a dipole magnetic field geome- 
try. Hence a determination of the magnetic field and its 
structure is critically important. 

Precise analysis of th e time domain of r adio emission 
from TVLM 513-46546 (iDoyle et alJlMtOft a nd sporadic 
radio emission from UCDs jAntonova et al.| [2007) show 
that large-scale fields may be stable on UCDs for long pe- 
riods, from a few months to years. The steady magnetic 
fields on UCDs are also confirmed by t he multi-frequenc y 
observations of a late-M dwarf binary (Os ten et al.ll2009l) . 
The field strength can be determined using two spe- 
cific radiation mechanisms - gyrosynchrotron or electron 
cyclotron maser (ECM). The first mechanism s uggests 
a fie l d stre ngth in the range of 0.1—1000 G (jBergerl 
120021 I2006D . while the latter implies a kG field. How- 
ever, the form of the frequency-field strength relation 
becomes complicated when the ECM mechanism is ap- 
plied to a many electron system since (1) the absorption 
and emission of different layers in the magnetosphere 
(or atmosphere) would be significant due to the differ- 
ent plasma environments and magnetic field configura- 
tion (e. g. see the discus s ion of the gyromagnetic absorp- 
tion in IMefrose fc Du lk (1982)); and (2) as we will show 
in this paper, for the motion of a group of electrons in 
a magnetic field, the ECM can generate a multiple peak 
structure for the spectral energy distribution. Hence, 
coverage of the full dynamic radio spectrum, including 
the low frequency band (hundreds of MHz) and the very 
high frequency band, are important for a proper under- 
standing of the radiation process. 

The observed power-law radio continuum and low level 
circular polarization (<40%) from several UCDs, such 
as for the M8.5 dwarf DENIS 1048 -3956 between 3 - 
30 GHz in four 2 GHz bandwidths (jRavi et al.|[20TTI ). 
may be interpreted as gyrosynchrotron radiation if the 
surrounding plasmas is optically thin. On the other 
hand, the high brightness temperature (~ 10 15 K) and 
highly (up to 100%) circul a r polarization of t he radio 
pulses (jHallinan et al.1 120071 : iBerger et all 120091 ) suggest 
that the dominant emission mechanism is the ECM. 
This mechanism was initially as sumed to be driven by 
a loss-cone velocity distribution (jMelrose fc Dulkl (1982) 
and references therein), but was subsequently devel- 
oped to ring shell distribution or horseshoe distribution 
(jPritchett fc Strangewavlll985| ). 

The operation of the ECM is rather simple, i.e. elec- 
trons with an anisotropic distribution transversely move 
in an external magnetic field. This leads to the ap- 
plication of the ECM to the radio emission from the 
solar pla n ets, m agnetic-chemically peculiar stars (e.g. 
lLo et al.l (|2012| )') to some comp a ct ex t ragala c tic ra - 
dio source s (e.g. IMeh^seJiDulk] (fl982l) : IPuIkl (fl985h ; 
ITreumannl (|2006l ) and references therein). The gener- 
ation of the auroral kilometric radiation on the Earth 
has been interpreted in terms of the ECM, where the 
velocity distribution of electrons may not be a loss- 
cone caused by the magnetic mirror effect; but due 
to a horseshoe distribution associated with the accel- 
erati on of particles in a magnetic field- aligned electric 
field (|Wu fc Led [1971: IChiu fc Schullfl^lErgun et al.1 
2000) . A similar interpretation can be applied to the de- 



camet ric radiation on Ju piter, Saturnian kilometric radi- 
ation (Zarka 1998, 2004) an d solar millisecond micr owave 
spikes (|Aschwandenlll990bl: iFleishman et al.ll2003D . The 
ECM can also be a strong candidate for the possible 
presence of radio emi s sion in exoplanets (iZarkal 120071 : 
iGriefimeier et al.1 120071 Uardine fc Cameron! 120081 ). Fur- 
thermore, it can be an effective mechanism for the radio- 
frequency heating of X- ray emitting plasma in solar flares 
(jMelrose fc Dulkl ll984[ ). Recently, it was suggested that 
the ECM generated by the low-density relativistic plas- 
mas in many fine localized regions can interpret the 
high brightness tempe rature detected from Blazar jets 
( Begcl man et al.ll2005D . More discussi on and app l icatio n 
of the ECM can be seen in a review in ITreumannl (|2006l) . 

In fact, gyrosynchrotron radiation and ECM belong 
to the same family - the motion of electrons in a mag- 
netic field. ECM may efficiently heat the surrounding 
electrons to form a high-energy tail or even a bump in 
the velocity space that induces gyrosynchrotron radia- 
tion to contribute to the radio continuum. Combin- 
ing the short time-scale and self- quenching features of 
the ECM, we can also understand the high brightness 
temperature of the radio pulses. Electron beams would 
be common in the context of the astrophysical process 
since there are plenty of sources for generating them. 
Recent cool atmospheric models indicate that collisions 
of significant volume of molecular clouds may trigger a 
tempestuous discharge process such as lightning, result- 
ing in a high-degree ionization in the local molecules or 
atoms, and the rel ease of a large number of electrons 
(|Helling et al.ll20Tl[l . which increases the probability of 
magnetic reconnection events. The electrons may be re- 
leased and accelerated from the magnetic reconnection 
or ou tflow jets indicated by oxygen forbidden emission 
lines (jWhelan et al.ll2007| ). which might result from the 
intense activity below the chromosphere of UCDs. 

In order to understand the radio emission from ultra- 
cool dwarfs and infer the magnetic field and the plasma 
environment, an investigation of a many electron system 
moving in an external magnetic field is essential. Numer- 
ical simulations can provide the opportunity to obtain 
the detailed process self-consistently and an interpreta- 
tion for the radio emission. In this paper, we attempt to 
interpret the radio pulses from UCDs using an electron- 
beam (or current-beam) driven ECM, with concentration 
on the microscopic energy transformation by treating the 
electron population as charged particles in a simulation 
box. We also investigate the growth rate and polariza- 
tion of the released electromagnetic (EM) waves and the 
spectral energy distribution (SED). In §[5j we briefly de- 
scribe the physical model, numerical method, and the 
initial conditions to carry out the simulations. In §[3] and 
§[U we present the results for the non-relativistic beam- 
driven and loss-cone-driven ECM. In §[5j we present the 
results for the relativistic beam-driven instability. We 
make a brief comparison with the observations in §[51 
then summarize the simulations and draw conclusions in 

§m 

2. CONFIGURATION OF SIMULATION 

2.1. Physical model 

We assume that electron beams are generated by some 
intense events on UCDs, e.g. magnetic reconnection or 
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jet events. When the electron beams move to the mag- 
netosphere of the UCDs, they interact with the magnetic 
held and the surrounding plasmas. In the present study, 
we neglect the influence of heavy ions. Here, we investi- 
gate the energy transfer, including the induced EM field 
energy, the drift kinetic energy and thermal kinetic en- 
ergy of electrons, plus the growth rate and polarization 
of the EM fields, and the spectral energy distribution. 

We start the simulations from the fundamental phys- 
ical laws. The EM fields and the interaction between 
them and electrons can be described by Maxwell's equa- 
tions, i.e. Ampere's law, Faraday's law of induction, 
Gaussian's law for magnetism, Gaussian's law and the 
definition of current 



V x B = p 3 



1 <9E 

~dt' 



V-B = 0, 



V-E= r. 



V- J = --f 



dp 



(i) 



(2) 



(3) 



(4) 



(5) 



where B and E are the magnetic field and electric field 
respectively, J the current, t the time, c the speed of 
light, p the charge density, eo the permittivity, po the 
permeability. 

The motion of electrons is governed by the Lorentz 
force, which can be written as 



q(E + v X B) = 



drrav 

~dT' 



(6) 



where v is the velocity of one individual particle, q the 
charge of a single particle (here it is for an electron), 
m the electron mass. In the non-relativistic case, we 
have m — m c where m c is the rest mass of the electron. 
In the relativistic case, we have m = ym e where 7 = 
— ; — i is the Lorentz factor. The relativistic case is 

also a general case for the motion equation of particles. 

2.2. Initial configurations of the simulations 

These equations were solved self-consistently as a pure 
initial value problem using a particle-in-cell method in 
a two dimensional space (x— and y~ direction) and 
three velocity and field dimensions (x—, y—, z— direc- 
tion). Some of the numerical meth ods used here ca n from 
lOmura fc Matsumotol (|1993l ) and lOmural (|2005l ). The 
Bunem an-Boris method was used to solve the equation of 
motio n (|Hocknev fc Eastwo od 1981; Bird sall fc Langdonl 
I1985T) . The equation of continuit y of charge was solved by 
a cha rge conservation method (jVillasenor fc Bunemanl 
[19921 ). The spacing grid for the electromagnetic field 
in the present simulations is 64x64. The time step 
in each simulation is At = 0.001 where v pc is 

the electron plasma frequency, while the space step is 
Ax = Ay = 0.125 Ad where Ad = Vth/(27rf po ) is the 
Debye length, with vth being the thermal velocity of the 



electrons. All the velocities in the simulations are nor- 
malized by the speed of light c. The charge-mass ratio of 
electron is assumed to be -1. We use an open boundary 
in the simulated system. This simulation configuration 
and the intrinsic properties of the electrons do not vary 
in any of the simulations. 

In each of the simulations, we assume a constant exter- 
nal magnetic field Bo exists in the spatial x — y— plane, 
and the angle between Bo and x— direction is defined 
as 6 with 0° ^ 9 ^ 90°. The charged particles are 
initially distributed in the x — y— plane randomly. We 
assume that background thermal electrons may exist in 
the radio emission region, i.e. the magnetosphere of ul- 
tracool dwarfs, with number density nth and thermal ve- 
locity Vth- The injected electrons have a number den- 
sity rid, thermal velocity v' th and drift velocity Vd along 
the x— direction. Figured] shows the spatial simulation 
box schematically. In the simulations, we determine the 
strength of the external magnetic field via its close re- 
lation with the cyclotron frequency v cc rj 2.8i?o MHz. 
The relation between plasma frequency and the number 
of electrons is u pc w 8.98 x 10 _3 (n o cm" 3 ) 1 / 2 MHz. 

Since we do not know the electron density in the radio 
emission region on ultracool dwarfs, a range of values and 
related plasma parameters are listed in Table [TJ From 
the values in the Table, we see that the present simula- 
tions are in a relatively micro-region and short timescale 
(0.1 ns - 10 fjs). 

In order to see the influence of the above parameters 
on the released EM waves, we set a group of standard 
values for them (see Table[2]). In this model, we assume 
the direction of the external magnetic field parallel to 
the y— direction. The cyclotron frequency is set to be 10 
times the plasma frequency. We take the thermal veloc- 
ity of the background electrons and the drift electrons as 
0.01c. This means that the temperature of the electrons 
is about 3xl0 5 K, determined by T = ^m e v^ h /k in the 
non-relativistic case, where k is the Boltzmann constant. 
We take I'd = 0.05c. We vary the value for one of the 
parameters whilst the other parameters remain as the 
standard values. These parameters and their values are 
summarized in Table[2] We will interpret these parame- 
ters in § 13.21 The standard values of these parameters are 
derived from estimations of solar bursts (typically 0.1c 
to 0.5c for the d rift ve l ocity and 0.002c to 0.05c for the 
thermal velocity, [Dulk (1985])), and the studies on auro- 
ral kilometric radiation on the Earth, Jovian millisecond 
bursts, and Saturnian kilometric radiation (~ 1 — 10 keV 
for the energe tic electrons a nd ~100 eV for th e ther- 
mal el ectrons. IZarkal (119981): IHess et all (I2007bl) ; IZarkal 
(|200l ; IHess et all (|2007alklLamv et al.l (^Oldf fh The val- 
ues of the electron veloc ities in the relativist ic case (see 
§[5]) refer to the work of lLouarn et al.l (|1986l ) and refer- 
ences therein. We choose the optional values over a wide 
range so that the approximate functions between energies 
and the parameters can be obtained. 

In this paper, we distinguish the irregular thermal mo- 
tion of the electrons and their uniform motion. The ther- 
mal energy of the electrons in the non-relativistic case is 
defined and calculated from the thermal motion of the 
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electrons by 



Eth — E t k — Ed 

n 1 

i=l 



1 



(7) 



where n is nth for the background electrons and n d for 
the drift electrons, m e the electron mass, Vj the velocity 
of the i th particle. On the right hand side of this equa- 
tion, the first term represents the total kinetic energy of 
the system E^, and the second term describes the drift 
energy of the electrons E d , 

In the relativistic case, the energy of one electron is 
defined by the energy-momentum relation 



Et 



P^ 2 

(7 m c 



2 4 



„2\2 



(8) 



where p io = 7m e Vi is the momentum of the i th electron, 
1 the Lorentz factor. Therefore the kinetic 
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Vl-(vi/c) 2 

energy of one electron can be expressed as 



E, 



ick 



E„ 



= (7- l)m c 2 . 
Then the total kinetic energy of the system E t ^ is 



E, 



tk 



ick- 



(9) 



(10) 



For the drift energy of the system Ed, we first define 
7d = „„ 1 , = . Then we have 



>/i-(£?=iV</nc)= 



E d = ( 7d - 1) 



(11) 



So the thermal energy is E t h = E t k — E&. The defini- 
tion of the drift energy of the system in both the non- 
relativistic and relativistic cases realises the uniform mo- 
tion of the system along one direction via eliminating the 
irregular random motion of the electrons. Under these 
definitions, the non-relativistic case (i.e. Eq.[7]) is a good 
approximation of the relativistic case (i.e. Eqs.[51 O flUI 
Hip when particles move with low velocities compared to 
the speed of light. Our calculations indicate that the 
equations for the relativistic case are valid for the ve- 
locity range Vi £ [0, c), while the non-relativistic case is 
only valid for low velocity particles. A general definition 
of temperature then becomes kT = (7 — l)m c c 2 . 

The energy of the EM field w can be described by the 
Poynting theorem, so that we have w — ^eoE 2 + 
where we take eo = 1 and fXo = The first term 

on the right hand side of this equation gives the electric 
field energy while the second term represents the mag- 
netic field energy. All kinds of energies are normalized 
by the initial total energy of the system to satisfy en- 
ergy conservation. We exclude the energy of the external 
magnetic field B . 

The growth rate is defined as 



2At 



(12) 



The degree of linear polarization in the EM fields II is 

w 1 — w u 



n 



(13) 



where w± is the EM field energy in the direction perpen- 
dicular to the external magnetic field, w u the EM field 
energy in the direction parallel to the external magnetic 
field. 

The background electrons are assumed to be in ther- 
mal equilibrium, so their velocity distribution obeys a 
Gaussian distribution. In this paper it is always taken 
as the following expression 



fob = n th exp(- 



-)• 



(14) 



1 Ui 



where v x , V y , v z are the components of the velocity Vb 



,y- 



direction 



of the background electrons along ir- 
respectively. 

The injected electrons have an intrinsic thermal ve- 
locity distribution but due to the acceleration they will 
obtain a drift velocity along some direction. For the pur- 
pose of simplicity, we assume all of the injected electrons 
are accelerated along one direction. We take the direc- 
tion as the x-direction in our simulations, so that the 
initial velocity distribution for the injected electrons can 
be expressed as 



fod = n d (- 



^reM- K ~ Vd)2 T y2+Vz ), (15) 



,'2 



"th 



1 th 



where / = 0,1,2,3,..., and v x , v y , v z are the compo- 
nents of the velocity v; n of the injected electrons along 
the x—,y—,z— direction respectively. I is the parame- 
ter which describes the size of the loss-cone in velocity 
space. In the present simulations, we take / = and 3. 
All the particles are assumed to be randomly distributed 
in space. 

3. RESULTS I: L=0 

When I = 0, Eq.[T5l becomes a Gaussian distribution. 
In this section, we investigate the influence of the in- 
jected electrons with the Gaussian velocity distribution 
on the evolution of the space and velocity distribution of 
the electrons, the energy conversion efficiency, the trans- 
fer between drift kinetic energy and thermal energy of 
the system, and the growth rate and polarization of the 
released EM waves. 

3.1. Standard model for beam-driven ECM 

We first present the simulation with standard parame- 
ters in detail; this we call the standard model. This stan- 
dard configuration means that the background electrons 
and the injected electrons have the same temperature 
(10 5 K) but the density of the latter is hig her. (Note we 
use the definition of the temperature in § 12.21 which re- 
flects the measure of the velocity dispersion) . This is for 
the purpose of modelling a dense electron beam injected 
into the magnetosphere of an ultracool dwarf. 

Figured] shows the spatial evolution of 
the electrons (A movie is available for the 
spat ial evolution of the system at URL: 
http: /7www.arm.ac.uk/~syu/2decm/beam/beam_space/[ ). 
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The time in each snapshot is t = 0, t = 0.3, t = 1.04, 
t = 2.86, t = 4.9, t = 6.52 from the top-left to 
bottom-right. Due to the existence of the external 
magnetic field, the motion of the electrons along the 
£— direction is confined and they can only freely move 
along the y— direction which is parallel to the magnetic 
field. The motion of the injected electrons as a whole 
should be in a helical orbit with the radius determined 
by = m °g d since the induced magnetic field is very 

small. In this simulation, is about 0.715. Since we 
only perform a 2D simulation, we see the oscillation of 
the injected electrons in the x — y-plane instead of a 
helical motion. 

The mix of background electrons and injected elec- 
trons in the velocity space triggered by the EM field is 
rather interesting. Fig. [3] shows snapshots of the velocity 
distribution of the electrons at the same time as in Fig. [5] 
(A movie is available for the velocity ev olution a t URL: 
http://www.arm. ac.uk/~syu/2decm/beam/beam_velocity/| ). 
As seen in the top-left panel in the figure, the velocities 
of both background electrons and injected electrons 
are initially a Gaussian distribution, with the injected 
electrons having a drift velocity of 0.05c. When time 
evolves, the current generated by the injected electrons 
induces a strong electric field along the direction per- 
pendicular to the external Bo which only alters the 
direction of the flow. This electric field accelerates a 
fraction of the background electrons so that a tail at 
high velocity is developed in v x and v z (see the velocity 
distribution in each time snapshot), whilst the injected 
electrons are decelerated by the electric field, losing 
their drift velocity along the a;— direction gradually and 
wavily. Also because of the perpendicular electric field, 
the injected electrons gain a drift velocity in v z which 
oscillates between —0.1c and 0.1c. The appearance of 
a double peak in the distribution of v y is due to the 
magnetic constraint and the acceleration of the induced 
electric field on the perpendicular velocity v x — v y . One 
can imagine that if the evolutionary time is sufficiently 
long the electrons may separate into two groups. One 
group having a velocity along the direction of the 
external B , while the other group has the opposite 
velocity direction. 

The consequence of this process is that the velocities of 
the electrons evolve from a concentrated Gaussian distri- 
bution to an expanded quasi-Gaussian distribution with 
a high velocity tail. During the diffusion process of the 
particles in velocity space, an energy transfer occurs be- 
tween the kinetic energy of the electrons and the induced 
EM field energy. As the injected electrons move in the 
external Bo at time t = 0, they start releasing their ki- 
netic energy to the EM wave energy in the manner of an 
increase in the induced EM field strength. After some 
time, the EM field energy will reach its maximum whilst 
Etk approaches a minimum. Then the EM field energy 
may be absorbed by the electrons to compensate for lost 
kinetic energy. The E t k will increase after a short time 
as the field energy decreases. Oscillations in E t k and the 
field energy will last for some time until the system is 
balanced and there is an anti-phase relation between the 
two kind of energies. 

Figure^] illustrates the evolution of the total kinetic 
energy E t k, drift energy Ed and thermal energy E t h of 



the electrons and the field energy w. We plot the same 
time points for the space and velocity distribution of the 
electrons in this figure as solid black circles. 

As we can see from this figure, E t k and w are exactly 
anti-phase as expected. The multi-peaks of w should 
be associated with the different wave mode in frequency 
space which will be shown later. The relation of the fine 
structure of Ed and E t h is not as obvious as the rela- 
tion between E t k and w since the interaction between 
E<i and E t h is via the EM field as a time and space de- 
lay can affect the phase relation. However, we see that 
-Ed decreases dramatically at the starting time when E t h 
increases rapidly. After about 1.5 J^, 1 , when Ed and 
E t h approximately have equal values, the decrease of Ed 
slows, and so does the increase in E t h- At time 7.1 v~ e }, 
at least 70% of Ed is converted to E t h- This is interesting 
because it means that the transverse motion of electrons 
in an external magnetic field may be an efficient way to 
heat the ambient electrons. 

The maximum growth rate of the EM wave in this 
simulation is about 9.68xl0 2 is pc . The polarization of 
the released EM waves is highly linear or circular, de- 
pending on the initial configuration. More details on the 
growth rate and polarization in the standard model will 
be shown in the next section. Figure[5] shows the evolu- 
tion history of the EM waves energy perpendicular and 
parallel to the external magnetic field. Again, we use the 
solid black circles to denote the time points for the space 
and velocity distribution of the electrons shown in Fig. [2] 
and [3] From this figure, we clearly see that most of the 
EM waves are polarized in the direction perpendicular 
to Bo, whilst only a very small fraction of the EM waves 
is released parallel to B (^10~ 3 of the perpendicular 
energy). 

3.2. Influences of parameters 

In this section, we investigate the influence of the fol- 
lowing parameters on the energy history, growth rate and 
polarization of the EM waves. 

(a) The strength of the external magnetic field Bq 
Bo can be determined by the cyclotron frequency 

B « 0.357 x 10 3 -^f- Gauss. (16) 
GHz 

To date, all detected radio emission of ultracool dwarfs 
are in the GHz band, while observations performed with 
the NRAO very large array in 2007 show no trace of 
radio emission from two UCDs at 325 MH z, placing an 
uppe r flux limit of ~ 900 ^Jy at 2.5<r level (paeger et al.l 
201 3). From this, we infer that the magnitude of Bq is 
few kilo- Gauss (if the radio emission is truly at the cy- 
clotron frequency). Here, we take f C e/^po = 0, 5 and 10 
to see how the magnetic field can affect the simulation 
results. When t'co/^'pc = 0, there is no external magnetic 
field. 

(b) The angle 9 between the magnetic field and the drift 
velocity 

6 is one of the crucial parameters to influence the di- 
rection of the radiated EM waves, i.e. the polarization. 
When 6 — 0°, the injected electrons move uniformly par- 
allel to the external Bo in addition to the irregular ther- 
mal motion. When 9 = 90°, the motion of the injected 
electrons is perpendicular to Bq. 



(> 



Yu et al. 



(c) The drift velocity Vd 

We expect that the radio emission may be enhanced 
by increasing the value of Vd- In this section, we take 
Vd = Oc, 0.005c and 0.01c to avoid the relativistic effect 
where gyrosychrotron emission plays an important role. 
When Vd = 0c, we see the effect of pure thermal electrons 
moving in the external Bo on the induced EM waves. 

(d) The temperature T 

We investigate the response of the radiated EM waves 
and the transfer between the energies by varying vth- 
We take «th = 0.005c, 0.01c, and 0.05c for which the cor- 
responding T is 7.5 x 10 4 , 3 x 10 5 , 7.5 x 10 6 K. 

(e) The background electrons by changing nth/^d 
Since our computation capability is limited, we only in- 
vestigate the effect of the existence of the background 
electrons on the EM waves and the energy transport. 
We take n t h/^d = 0,1,2. When n t h/«d = 0, there are 
no background electrons. 

Figure[5] shows the energy evolution in each simulation. 
We see that the EM field energy (left panels) are in a 
very low level in all three specific cases (see the colored 
lines) . In the first case, there is no magnetic field (which 
can be understood as the injected electrons just pass by 
very quickly without sufficient energy transfer via the 
EM field). In the second case, the motion of the injected 
electrons is parallel to the external magnetic field. In 
this case, the induced electrons can not remain within the 
simulation box since the magnetic field can not constrain 
them. In the third case, the drift velocity of the injected 
electrons is 0c; although the injected electrons can remain 
and interact with the background, there is no coherent 
current, i.e. the field energy is still small although higher 
than in the other two cases. 

A common expression of these cases is v x Bo = which 
is the Lorentz force induced by the external magnetic 
field. This force makes the initially coherent current bend 
(i.e. the electrons with drift velocity), leading to spatial 
curled EM fields that is the medium to accomplish energy 
transport among different kind of energies. 

In other cases except the above three cases, the field 
energy can maintain a much higher level after they reach 
the maximum, typically orders of 2 — 3 that of the above 
cases. As shown in Fig. [5J the field energies in all cases 
oscillate with large amplitude caused by the transfer be- 
tween the kinetic energy and field energy. In other words, 
these oscillations reflect the emission and absorption of 
electrons to the EM waves. With the achievement of the 
diffusion process of the electrons in the velocity space, a 
dynamic balance is approached. This results in a gradual 
decrease in the amplitudes of the oscillations with time. 
The time for relaxation of the field energy is >10 
which is much longer than the time taken for the field 
energy to reach maximum, ^0.14 v~£. 

As expected, increasing Bq by a factor of two, i.e. 
Vce/vpe from 5 to 10, leads a rise in the field energy. 
It seems that the mean value of the field energy is not 
sensitive to 9 when 9 = 45° and 90°. However, the 
modes (or direction) of the EM waves are affected. In 
the case of 9 = 45°, the EM waves have a similar en- 
ergy level in the direction perpendicular and parallel to 
Bo, which is shown in Fig. [7] where we discuss polariza- 
tion. The increase of Vd from 0.05c to 0.1c only rises the 
energy level slightly in the present simulations. In fact, 



when Vd is sufficiently high, we have to consider the rela- 
tivistic effect and hence gyrosynchrotron radiation which 
will be addressed in §[5j The mean energy level of the 
EM waves is not sensitive to the thermal velocity and 
background electron number density. Note that in the 
standard model, the initial density of the background 
electrons is ~ 100 times less than that of the injected 
electrons. 

The polarization of the EM waves, or their energy dis- 
tribution with respect to the magnetic field direction, 
is shown in Fig. [7] quantitatively. The EM waves in 
the standard model are 100% linearly polarized. When 
v x Bo = 0, the EM waves frequently switch their direc- 
tion from parallel-dominant to perpendicular-dominant, 
and rarely do they have up to 50% linear polarization. 
This behaviour indicates that the waves have a significant 
linearly-polarized component. In addition to 9, thermal 
motion (i.e. temperature) of the electrons can influence 
the level of polarization. As we see from the bottom 
second panel in Fig. [3 the electrons with temperature 
7.5 x 10 6 K can generate the EM waves with ~50% — 
^90% linear polarization, while 100% linearly polarized 
waves are generated by the electrons with temperature 
7.5 x 10 4 and 3 x 10 5 K. A low density of background 
electrons does not alter the highly linear polarization in 
the present simulations. 

The dissipation of the drift energy in the injected elec- 
trons (when they move in the external magnetic field) 
is important because this may be sufficient to increase 
the thermal energy of the system. We show the history 
of the drift energy and thermal energy in Fig. [5] Com- 
paring the right panels in Fig.[BJ we find that: (i) when 
v x Bo = 0, the energy transfer is the least efficient. In 
this case, there is almost no energy exchange; (ii) when 
initially vxBo ^ and Ed > E t h, Ed can be transported 
to -Eth rapidly; the timescale in which Ed and Eth reach 
the same value is about ~ 1 — 5 v~} . 

Figure|8] illustrates the growth rate of the EM waves in 
the simulations. The increase of the system temperature 
(3 x 10 5 — 7.5 x 10 6 K) will suppress the growth rate 
significantly. The number of the electrons also affects 
the growth rate, but their relation is not so obvious. 

3.3. Spectrum 

The EM field (wave) energy is the integral (or sum if 
the signal is discrete) of the contribution from different 
frequencies. In order to obtain the energy distribution of 
the EM field in frequency space (i.e. the spectral energy 
distribution, SED), we perform a fast Fourier transform 
(FFT) to the EM field energy history. Figure[9] illustrates 
the SED. 

We clearly see many emission and absorption lines in 
Fig-El which represent the EM field energy distribution 
at different frequencies. We find the majority of the field 
energy is from frequencies < 80^ pc with bandwidth at 
half maxima ~ 6^ pc . 

4. RESULTS II: L = 3 

In order to see the interactions of the electrons with 
different initial velocity distributions, we perform an- 
other series of simulations with 1 — 3. In this case, 
when Vd — 0, Eq.[l5] is a typical loss-cone velocity dis- 
tribution. Note that in this section, we take Vd = 
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in the standard model to exclude the effect of the 
electron-beam. We vary both the thermal velocity 
of the electrons to see the effect of the temperature 
on the energy exchange and release (Movies are avail- 
able for t he spatial evolut ion of th e system at URL: 
http:/ /www. arm. ac.uk/~syu/2decm/losscone/ 
lossco he~l=3-space/ and for the velocity evolu tion URL: 
http://www.arm.ac.uk/~syu/2decm/losscone/ 
lossconeJ=3_veiocity/j. 

4.1. Standard model and influence of parameters in 
loss-cone-driven ECM 

Figurc[lO] illustrates the EM field in the left panels and 
in the right panel we have the thermal and drift energy 
of the electrons from the loss-cone driven ECM. From 
this figure, we see that only a very small fraction of ki- 
netic energy is converted to EM field energy if there is 
no external magnetic field or the loss-cone velocity dis- 
tribution is along the direction parallel to the external 
magnetic field. 

In other cases, an increase in any one of the param- 
eters Bo, 9 and v' th leads to an increase in the induced 
EM field energy. Rising the temperature of the back- 
ground electrons (i.e. the thermal level) can suppress 
the growth of the field energy. Since we only vary the 
number of background electrons in a very small range, it 
does not affect the field energy significantly. As seen from 
the left-bottom panel in Fig.[K)l increasing the number of 
background electrons decreases the induced EM field en- 
ergy. These results are consistent with the growth rate of 
the EM wave illustrated in Fig.[TTJ From the right pan- 
els in Fig.HUl we see that the drift energies are rapidly 
dissipated, eventually leading to irregular motion of the 
electrons in the external magnetic field. 

The history of the degree of the polarization of the EM 
waves in each simulation is shown in Fig. [12] The condi- 
tions for the circularly polarized EM waves are notable, 
ranging from (1) decreasing the magnetic field; (2) vary- 
ing the angle 9 from perpendicular to non-perpendicular; 
to (3) increasing the thermal level of the background elec- 
trons. We suggest these conditions may be associated 
with the circularly polarized components of the radio 
emission from ultracool dwarfs. 

4.2. Spectrum 

The influence of the parameters on the spectrum is 
shown in Fig. [13] in which we find that the magnetic field 
and the angel 9 can affect the spectral energy distribution 
significantly, while the thermal effect and the number 
of background electrons only plays a minor role. Some 
frequency bands are notable, e.g. from 10 to 70f pc . 

4.3. Comparison between beam-driven and 
loss-cone-driven ECM 

Comparison between the different initial velocity dis- 
tributions can help us understand the roles of the initial 
parameters in the process of releasing EM wave energy 
and the transfer between drift energy and thermal en- 
ergy. Combining Fig. |6]& [7] (beam-driven) and Fig.lTOlfe 
[T2l (loss-cone-driven), we see that: 

1) the existence of the drift kinetic energy of the elec- 
trons can be considered as a coherent current (which is 
necessary to generate the intense EM field energy) , while 



the form of the initial velocity distribution is not impor- 
tant; 

2) in order to efficiently obtain the intense EM field 
energy, the external magnetic field plays a crucial role. 
The angle between the magnetic field and the coherent 
current significantly affects the strength of the released 
EM field energy in a non-linear relation (see Fig. [5] and 
[TT]) and also the propagation direction of the EM waves; 

3) pure thermal motion of the injected electrons in the 
external magnetic field may play a role in generating the 
EM waves (e.g. when Vd = in the beam-driven ECM, 
middle panel in Fig. [6]), while the thermal level of the 
background electrons mainly suppress the generation of 
the EM waves. 

The SED of the beam-driven and loss-cone-driven 
ECM indicate that all the parameters (except the num- 
ber of background electrons in the present simulations) 
can affect the SED, however certain harmonic frequency 
bands will appear if the coherent current is sufficiently 
strong, for example see the region from 10 to 70^ pc . 
There is a negligible signal in the very high frequency 
band > 100^ pc . Also, it seems that the SED weakly de- 
pends on the number of background electrons, but this 
needs further investigation since we do not vary the num- 
ber of particles over a sufficiently wide range in the sim- 
ulations. 

5. RELATIVISTIC BEAM-DRIVEN INSTABILITY 

In this section, we show the case where a relativistic 
electron-beam moves in an external magnetic field. We 
take I = in Eq.[15]and set the same standard parameter 
values as in §[3] except here we take a larger value for the 
drift velocity of the electrons, i.e. Vd = 0.98c for the 
standard model; 0.58c and 0.78c for the optional values. 
Furthermore, we only vary the thermal velocity of the 
background electrons, with values of Vth — 0.001c, 0.01c 
(standard value), and 0.1c. 

5.1. Influences of parameters 

The history of the three kind of energies, the polar- 
ization of the induced EM field energy and their growth 
rates are illustrated in Fig.Q3] [T5] and [THJ It is seen that 
the background electrons in the simulations have a neg- 
ligible influence on the induced EM field energy and the 
linear polarization while the thermal energy level of the 
background electrons can suppress the growth of the EM 
waves. 

From the left panels in Fig.Q3] we find that the direc- 
tion and values of the magnetic field and the drift velocity 
can significantly affect the induced EM field. The exis- 
tence of the magnetic field is necessary in order to obtain 
the fast growth of the EM field energy and the energy 
conversion efficiency seems to have a non-linear relation 
with the magnetic field strength. 

The right panels in Fig.[T4l illustrates the drift energy 
and the thermal energy of the electrons. If the magnetic 
field strength is sufficiently high to constrain the motion 
(or spatial position) of the injected electrons, the drift 
energy can be rapidly converted to the thermal energy in 
a time scale which is similar to the non-relativistic case. 
If the magnetic field is not very strong, we see that it is 
possible for the system to still retain some residual drift 
energy when the electrons escape from the local region. 
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The influence of the parameters on the polarization in 
the relativistic case (Fig. [T5|) is different from the non- 
relativistic case. The strength of the magnetic field in 
the present simulations affects the degree of the linear 
polarization. The existence of a strong magnetic field is 
essential to generate the linear polarization. The non- 
perpendicular angle between the magnetic field and the 
drift velocity plays an important role in the generation of 
the linear polarization. In fact, at high frequency, linear 
polarization becomes dominant in the relativistic case. 
The drift velocity, the thermal level of the background 
electrons and the relative number of drift electrons and 
background electrons do not affect the polarization sig- 
nificantly in the present simulations. As shown in Fig. ll6[ 
the growth rate of the EM field rapidly increases with the 
magnetic field strength and injection angle. However, it 
seems that in the relativistic case the increase of the drift 
velocity slightly decreases the growth of the EM field per- 
haps because of a strong interaction between the induced 
EM field and the high energy electron beam. The ther- 
mal level of the background electrons can suppress the 
growth of the EM field, and there is an ambiguous rela- 
tion between the number of background electrons and the 
growth rate of the EM field in the present simulations. 

The influence of the direction between the magnetic 
field and the drift velocity is very interesting since the 
maximum energy conversion efficiency occurs in the 
non-perpendicular injection which differs from the non- 
relativistic case. In order to determine the direction 
where we can obtain the maximum energy conversion effi- 
ciency, we have done some extra computations by varying 
the direction of the magnetic field. Figure[T7] illustrates 
the results in which we see that the maximum energy 
conversion efficiency occurs at ~ 75°. 

An interesting phenomenon in our simulations is that 
the velocity distribution of the beam electrons may 
evolve from an initially drifted Gaussian (or Gaussian- 
like) distribution to ring distribution (or incomplete- 
ring or spiral-ring or horse-shoe, e.g. see the standard 
models in both the relativistic case and non-relativistic 
case.) (A movie is available for the velocity evolution 
of the electron beam in the re lativistic case at URL: 
http : / / www, arm, ac . uk/~syu / 2decm / relativistic_beam / 1 
relativistic_beam_velocity/) or spherical-shell distribu- 
tion (or incomplete shell, e.g. if we change the angle 
between the magnetic field and the beam electron from 
perpendicular, e.g. 90°, to non-perpendicular, e.g. 45° 
or 75°.) This may imply that the ring and shell (or 
ring-like and shell-like) velocity distribution would have 
a common origin — beam distribution. The influence of 
the angle on the evolution of the velocity of electrons 
is mainly caused by the external magnetic field and the 
self-induced EM fields. The spatial current induced by 
the motion of the charged particles plays an important 
role in the process, which initially causes the variation 
of a spatial magnetic field and time-dependent electric 
field, i.e. Eq.[TJ This electric field will accelerate or de- 
celerate the electrons to deform the distribution values 
of the electron velocity. 

5.2. Spectrum 

The SED of the radiation in the relativistic electron 
beam apparently differs from that in the non-relativistic 
case. FigurefTBl illustrates the SEDs in the relativistic 



case. In the standard model, we clearly see that the en- 
ergy may distribute not only in the range of 10 to 100^ pc , 
but extend up to a frequency of ~ 500z/ pe and some neg- 
ligible signal at an even higher frequency harmonic. 

For a single relativistic electron, the peak value of 
its synchrotron radiation m ay be at v pea ^ 7 2 z/ ce sin# 
(|Rvbicki fc Lightmanl 11979] ). For a many electron sys- 
tem, the frequencies where we can obtain emission lines 
is strongly associated with the distribution of the Lorentz 
factor 7. FigurelT9l illustrates the distribution of 7 in dif- 
ferent parameters. In the standard model, we find that 
when the instability reaches saturation the distribution 
of 7 expands only a little in the high energy part. How- 
ever, when varying the angle 9 from 90° to 75°, we see a 
distinct high energy electron tail around 7 = 9.5 which 
corresponds to the kinetic energy = (7 — l)m c 2 « 
4.35 MeV. These high-energy electrons may contribute to 
possible X-ray emission from UCDs via thermal or non- 
thermal bremsstrahlung or even inverse Compton scat- 
tering. 

6. COMPARISON WITH OBSERVATIONS AND 
OTHER STUDIES 

It is possible to approximately compare our simulation 
results with the observed radio spectrum from UCDs 
using some reference values of v pc in Table[TJ For ex- 
ample, if we assume that the observed radio emission 
from TVLM 513-46546 at 8.5 GHz corresponds to the 
simulated peak frequency at ~ 12i^ pc in the case of 
^'ce/t'pc = 10 (the strongest in the simulation) in the 
non-relativistic beam-driven instability, then v po — 0.708 
GHz and v cc — 7.08 GHz. Thus the corresponding lo- 
cal plasma density is ~ 6.2 x 10 9 cm -3 and the strength 
of the magnetic field is ~ 2529 Gauss. We may find 
some weaker signal at a high frequency harmonic, such 
as - 30zv = 21.24 GHz, - A2v pc = 29.74 GHz, 
~ 50i/ po = 35.4 GHz. 

In the relativistic case, the first significant signal ap- 
pears at ~ 20^p G when v cc jv pc = 10, which might cor- 
respond to the observed radio signal at 8.5 GHz. This 
gives v pc = 0.425 GHz and v cc = 4.25 GHz. Thus the 
corresponding local plasma density is ~ 2.2 x 10 9 cm -3 
and the strength of the magnetic field is ~ 1518 Gauss. 
We may find some weaker signal at a higher frequency 
harmonic, such as ~ 40^ pc = 17 GHz, ~ 118^ pc = 50.15 
GHz, ~ 195fp C = 82.875 GHz, etc. This means the spec- 
trum may extend to the extreme high frequency band or 
even far infrared in the extremely relativistic case. 

If we choose v cc jv pc = 5 in the relativistic case, and 
assume that the first significant signal appearing at ~ 
llO^po corresponds to the observed radio signal at 8.5 
GHz, it gives u pe = 0.0773 GHz and v ce = 0.386 GHz. 
Thus the corresponding local plasma density is ~ 7.4 x 
10 7 cm~ 3 and the strength of the magnetic field is ~ 138 
Gauss. 

The above estimations of the properties of the emission 
region and the spectrum distributions in the presen t sim- 
ulatio ns ar e consis t ent w ith the analytical results in lDulkl 
<fl985h and lGlidel (|200l . 

We do not draw conclusions on the values of plasma 
parameters in this paper because firstly many parame- 
ters are free in our simulations. Especially two of the 
key parameters, the local electron density and magnetic 
field, are ambiguous. Secondly the present simulations 
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are confined in a local micro-region with no consideration 
of absorption and emission by other regions. The absorp- 
tion and re-emission by the plasma in other regions needs 
to be further investigated. The radio spectrum can be 
affected by the large-sc ale structure of a magnetic field 
(jKuznetsov et al.l 120 111 ) . A combination of observations 
at other wavelengths will help determine the configura- 
tion of the magnetic field and the plasma environment, 
and thus entirely understand the magnetosphere and at- 
mosphere of ultracool dwarfs. 

The simulation results of the radio pulse s from one 
of the radio active UCDs in lYu et all (120111) show that 
the loss-cone-driven ECM in lAschwandenl (| 1990a! ) can 
result in the release of ~ 0.5% of the kinetic energy, which 
therefore can generate a strong radio flux of up to a few 
mjy. This depends on other parameters, for instance the 
local plasma density. By comparison, the non-relativistic 
beam-driven ECM in the present model can release ~2% 
of the kinetic energy when the process reaches saturation, 
which is about 4 tim e s high er than the results in the 
model of lAschwan dcn (19903). 

Our simulation results are also comparable with 
the studies fo r auro ral kilometric radiation regions. 
IPritchett et al.l (|2002t ) investigated the auroral kilometric 
radiation source cavity using two-dimensional particle-in- 
cell simulations. They found the energy conversion effi- 
ciency is ~ 1 — 5%, depending on the velocity of the hot 
electrons (~ 1 — 10 keV). The time scale of the radiation 
bursts is ~0.5 ms, which is pro bably associated with th e 
local plasma density ( <1 cm~ 3 .IStrangewav et al.l (fl998h : 
iPerraut et all (|1990l ) ; fCalvertl dl98lD) and the weak mag- 
netic field. The energy conversion efficiency in our re- 
sults are in line with those in their results in the non- 
relativistic cases. The spectra obtained in the present 
simulations are also in good agreement with the previ- 
ous resul ts (IMelrose fe Dulkl 119821: lAschwandenl Il990at 
IPritchett! 119861 11984) - the waves are produced mainly 
near the electron cyclotron frequency. 

However, certain differences exist between our simula- 
tion results and others. Firstly, in our simulations, we 
can see the different oscillation modes from the energy 
time history which reflect the interaction between the 
particles and the induced electric field. Secondly, we 
find that the energy conversion efficiency can be influ- 
enced significantly by the drift velocity of the particles 
and the angle between the drift velocity and the mag- 
netic field, which is perhaps caused by the self-induced 
electromagnetic damping. The initial conditions and the 
technique to solve the Maxwell's equations in our sim- 
ulations and others are also different. We considered a 
spatially localized electron population with a beam or 
beam- like velocity distribution. In contrast, previous 
works c onsidered ring-like, horseshoe, and DGH d istri- 
butions (|Pritchett et al.ll2002t IPritchett! fl986l I1984D and 
a larger space structure for the magnetic field. The time 
advancement of Maxwell's equ a tions was performed in 
Fourier space in IPritchett et al.l (|2002l ). Instead, we used 
a staggered leap-frog method to solve the equations in 
the time domain with the positions and velocities of the 
electrons generated by a Monte Carlo method. 

Recent comput er simulation s and expe rimental lab- 
orator y work by I Cairns et al.l ()201lD and iVorgul et al.l 
(|2011D show that in the non-relativistic (or weakly rel- 
ativistic) case the energy conversion efficiency is 1—2%, 



which is consistent with our simulations. However, this 
energy conversion efficiency varies significantly in our rel- 
ativistic case, depending on the parameters, e.g. the 
angle, the drift velocity, the strength of external mag- 
netic field. It may be mainly caused by the effect of 
synchrotron radiation. 

7. SUMMARY AND CONCLUSIONS 

In this paper, we present the numerical simulations 
for electron-beam-driven and loss-cone-driven ECM with 
different plasma parameters and different magnetic field 
strengths. We find that the beam-driven ECM can be an 
effective mechanism to release EM waves and heat the 
surrounding plasmas. From the diffusion process of the 
electrons in velocity space and the energy distribution 
(see Fig. [3] and Fig. [19]) , a high-energy tail of the elec- 
trons may be rapidly developed along the direction near 
perpendicular to the magnetic field, which can eventu- 
ally evolve to moderately or strongly relativistic electrons 
depending on the initial energy of the electron current, 
and contribute towards gyrosynchrotron or synchrotron 
radiation. This may lead to the appearance of a radio 
continuum and the deformation of the SED. Also, these 
high-energy electrons may be important to generate X- 
ray emission. 

The computation of the degree of polarization indicates 
that the thermal level of the electrons can significantly 
affect the degree of the circular or linear components of 
the observed radio waves. In the case of the beam-driven 
ECM, the angle between the direction of the magnetic 
field and the injection direction of the injected electrons 
is another crucial factor to affect the degree of circular 
polarization in the radio waves. 

The SEDs of the radio waves depends weakly on the 
form of the velocity distribution of the electrons in the 
present simulations. The existence of the external mag- 
netic field and the angle between the direction of the 
magnetic field and the moving direction of the electron 
current can significantly affect the SED. Certain fre- 
quency bands, e.g. 10 to 70^ pc in the non-relativistic 
case and 10 to 600v pe in the relativistic case may appear, 
which increases the difficulty of finding the fundamental 
cyclotron frequency in the observed radio frequencies. It 
is however possible that magnetic field inhomogeneities 
may smooth out some of these bands thus producing a 
continuous spectrum. In order to determine the plasma 
frequency and the cyclotron frequency, wide frequency 
bands should be covered by future radio observations. 

The present study is limited in that only two- 
dimensional electromagnetic simulations are performed 
with no consideration of the change of the magnetic field 
configuration and the influence of gravity. We will con- 
tinue to develop the simulations to match the plasma 
environment and the magnetic topology on UCDs in or- 
der to understand their radio emission. 
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Drift electrons 




Figure 1. The schematic diagram to present the simulation box for the radio emission from ultracool dwarfs. In an orthogonal reference 
system with spatial direction x, y and z, ti x , v y and v z represent the velocity components of particles, while E x , E y , E z , B x , B y and B z 
represent the electric and magnetic fields components. 




Figure 2. The history of the distribution of the background electrons (red points) and beam electrons (blue points) in spatial x-y-plane. 
Different panels from top-left to bottom-right represent the snapshots for time t = 0.00, 0.30, 1.04, 2.86, 4.90 and 6.52 Up e . 
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Figure 3. The history of the velocity distribution of the background electrons (red points) and beam electrons (blue points). Different 
panels from top-left to bottom-right represent the snapshots for time t = 0.00, 0.30, 1.04, 2.86, 4.90 and 6.52 Vp e . 
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Figure 4. The energy history of the standard model. Blue line in the left panel denotes the total kinetic energy, and red line stands for 
the field energy. In the right panel, purple line is for the total thermal energy of the system, and the green line is for the drift energy of 
the injected electrons. Black points denote the time as in Fig.[2]and[3] 
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Figure 5. The history of the field energy along the direction of perpendicular (red line) and parallel (blue line) to the external magnetic 
field. Black points denote the time as in Fig,l2landl3l 
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Figure 12. The time evolution of the degree of linear polarization in the loss-cone-driven ECM. Each panel is for one parameter, i.e. from 
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Figure 14. Same as Fig. [6] but for the relativistic beam-driven instability. 
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Figure 15. The evolution of the degree of linear polarization with time in the relativistic beam-driven instability. Each panel is for one 
parameter, i.e. from top to bottom, v ce /v pe , 8, ly, v t ^, n^/n^. Different colors in each panel are for different values of the parameters as 
shown in the panel. 
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Figure 17. The energy history for the standard model with different angle between the external magnetic field and the injected velocity 
direction in the relativistic beam-driven instability, i.e. 45°, 55°, 65°, 75°, 85°, 90°. 
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Figure 18. Same as Fig. [9] but for the relativistic beam-driven instability. 
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Figure 19. The energy distribution of the electrons in the relativistic beam-driven instability. Top panel is for the initially mildly 
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Table 1 

Reference values of electron densities and related parameters. 
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Table 2 

Various parameters in the simulations (§03 and their values. Note that part of the values may vary in §[4]and §[5] 



Parameters 


Symbols 


Standard 


Optional 






values 


values 


Plasma trequency 




1 




Cyclotron frequency 


fee 


10 


0, 5 


Angle: Bo and 





90° 


0, 45° 


Thermal velocity 








of background e 


Vth 


0.01c 


0.005c, 0.05c 


Thermal velocity 








of drift e 


«£h 


0.01c 


0.005c, 0.05c 


Drift velocity 


v a 


0.05c 


0.0c, 0.1c 


Number of superparticles 








of background e 


"th 


16384 


0, 32768 


Number of superparticles 








of drift e 


n d 


16384 





